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ABSTRACT 

A multifrequency campaign on the BL Lac object PG 1553+113 was organized by the 
Whole Earth Blazar Telescope (WEBT) in 2013 April-August, involving 19 optical, 
two near-IR, and three radio telescopes. The aim was to study the source behaviour 
at low energies during and around the high-energy observations by the Major At¬ 
mospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes in April-July. We 
also analyse the UV and X-ray data acquired by the Swift and XMM-Newton satel¬ 
lites in the same period. The WEBT and satellite observations allow us to detail the 
synchrotron emission bump in the source spectral energy distribution (SED). In the 
optical we found a general bluer-when-brighter trend. The X-ray spectrum remained 
stable during 2013, but a comparison with previous observations suggests that it be¬ 
comes harder when the X-ray flux increases. The long XMM-Newton exposure reveals 
a curved X-ray spectrum. In the SED, the XMM-Newton data show a hard near-UV 
spectrum, while Swift data display a softer shape that is confirmed by previous HST- 
COS and lUE observations. Polynomial fits to the optical-X-ray SED show that the 
synchrotron peak likely lies in the 4-30 eV energy range, with a general shift towards 
higher frequencies for increasing X-ray brightness. However, the UV and X-ray spec¬ 
tra do not connect smoothly. Possible interpretations include: i) orientation effects, 
ii) additional absorption, iii) multiple emission components, and iv) a peculiar energy 
distribution of relativistic electrons. We discuss the first possibility in terms of an 
inhomogeneous helical jet model. 

Key words: galaxies: active - BL Lacertae objects: general - BL Lacertae objects: 
individual: PG 1553+113 
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1 INTRODUCTION 


BL Lac objects, together with flat-spectrum radio quasars, 
form the “blazar” class of active galactic nuclei (AGNs). 
Their observed properties, such as strong variability at all 
frequencies, high and variable polarization, apparent super¬ 
luminal motion of the radio components, are explained as 
due to beamed emission from a relativistic jet closely a ligned 
with the line of sight (e.g. lUrrv fc PadovaiiH ll995l L The 
broad-band spectral energy distribution (SED) of these ob¬ 
jects shows two bumps: the low-frequency one (from radio 
to optical/X-rays) is ascribed to synchrotron emission, while 
the bump at high energies (X-rays to 7 -rays) is likely pro¬ 
duced by inverse-Compton scattering of soft photons off the 
same rela tivistic electr ons responsible for the synchrotron 
radiation (lKonig]lll98M 

The source PG 1553-1-113 is one of the so-called “high- 
energy-peaked” BL Lacs (HBL), whose synchrotron peak in 
the SED typically falls in the UV-X-ray frequency range 
ijPalomo fc Trev3ll99fj l. The redshift of PG 1553-1-113 is 
still unknown. Many authors tried to set lower and up¬ 
per limits in an indi rect way, from the lack of the h ost 
galaxy detection (e.g. iTreves. Falomo fc Usle nghill2007l l or 
from the 7 -ray spectru m (e.g. lAbdo et aLlbOld L In contrast, 
IPanforth et al.l (| 2 fllfll j derived Arm constraints 0.43 < 2 < 
0.58 from the direct analysis of inter-galactic absorption fea¬ 
tures visible in the UV spectra acquired by the Hubble Space 
Telescope/Cosmic Origins Spectrograph (HST-COS). 

Very high energies (VHE, E > 100 GeV) observations 
with MAGIC in 2005-2009 showed modest flux variabil- 
ity, p ossibly correlated with the optical one (|Aleksic et al.l 
[2013). Ill 2012 the source was found in a flaring state at 
VHE and X-rays, b ut not at the GeV energies covered by 
the Fermi satellite (|Aleksic et al.l[2015l L New MAGIC ob¬ 
servations were performed in April, June and July 2013 
during a low emission state. This time, a broad multiwave¬ 
length campaign was organized by the Whole Earth Blazar 
TelescopeEl (WEBT), an international collabor ation born in 


1997 to study blazar em i ssion variability (e .g. Villnta^_et_^ 
2 OO 2 I: iRaiteri et al.l 200 (il: IVillata et al. 200(il:lB5ttcher et al 


2007i : [Larionov et al. 1 200ft R.aiteri et al.l 200911 . In this paper 


we present the results of the WEBT monitoring campaign, 
which was extended to all the 2013 observing season, and 
explore the synchrotron part of the source spectral energy 
distribution (SED). To this aim, we complement the radio- 
to-optical WEBT observations with the UV and X-ray data 
acquired by the Swift and XMM-Newton satellites in the 
same period. The results of the MAGIC observations will 
be analysed in a forthcoming paper (Ahnen et al. 2015, in 
preparation). 

The paper is organized as follows. In Sect. 2 we present 
the WEBT light curves in the optical BVRI and near-IR 
JHK bands as well as radio light curves at three wave¬ 
lengths. We also derive calibration of an optical photometric 
sequence in the BVRI bands and perform colour analysis. 
The Swift and XMM-Newton observations are analysed in 
Sect. 3 and 4, respectively, while in Sect. 5 the SED shape 
and variability is discussed with the help of past data by the 


^ Al ternatively, a hadron ic cascade scenario has been proposed 
fe.g.. [Bdttcher et alj|20l3 . and references therein) 

^ http://www.oato.inaf.it/blazars/webt 


Table 1. The optical, near-IR and radio observatories participat¬ 
ing in the WEBT campaign. 


Observatory 

Country 

Bands 

Abastumani 

Optical 

Georgia 

R 

Belogradchik 

Bulgaria 

BVRI 

AstroCamp 

Spain 

R 

Crimean 

Russia 

BVRI 

Michael Adrian 

Germany 

BVRI 

Mt. Maidanak 

Uzbekistan 

BVRI 

New Mexico Skies 

USA 

R 

Plana 

Bulgaria 

BVRI 

Rozhen^ 

Bulgaria 

BVRI 

San Pedro Martir 

Mexico 

R 

Siding Spring 

Australia 

R 

Skinakas 

Greece 

BVRI 

St. Petersburg 

Russia 

BVRI 

Teide 

Spain 

R 

Tijarafe 

Spain 

R 

Valle d’Aosta 

Italy 

BVRI 

ASV2 

Serbia 

BVRI 

Campo Imperatore 

Near-infrared 

Italy 

JHK 

Teide 

Spain 

JHK 

Medicina 

Radio 

Italy 

8 GHz 

Metsahovi 

Finland 

37 GHz 

Noto 

Italy 

43 GHz 


^ Three telescopes 
^ Astronomical Station Vidojevica 


Swift, the International Ultraviolet Explorer (lUE) and the 
HST satellites. Conclusions are drawn in Sect. 6 . 


2 WEBT DATA 

Optical observations for the WEBT campaign were per¬ 
formed with 19 telescopes in 17 observatories around the 
globe. These are listed in Table [1] with indication of their 
observing bands. 

Calibration of the source magnitude was obtained with 
respect to the reference stars 1, 2, 3 and 4 shown in Fig. 
[ 1 ] Photometry of Stars 3 and 4 was derived from the Sloan 
Digital Sky Su rve^B (SPSS), usin g uqr iz to UBVRI trans¬ 
formations by IChonis fc Gaskelll (I2OO8II . The SPSS magni¬ 
tudes for Stars 1 and 2 are affected by saturation, so we 
calibrated these stars through differential photometry with 
respect to Stars 3 and 4 using our best-quality data. The re¬ 
sults of our photometry are re ported in TablelH They a re in 
good agreement with those of lAndruchow et al.l (|201ll '> and 
jPoroshenko et al] ll2005l . [2014h . while the difference with the 
values contained in the USNO-A2.0 catalogue is up to sev¬ 
eral tenths of mag. 

In total, we collected 3908 optical data points in the 
period 2013 January 25 to September 11. Light curves in 
BVRI bands were built by carefully assembling the data sets 
coming from the different telescopes. Cleaned light curves 


® http://www.sdss.org/ 
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Table 2. Magnitudes (and errors) of the optical photometric sequence shown 
in Fig. [T] for calibration of the PG 1553+113 magnitude. The values de- 
rived from ground-based obse rvations by the WEBT for the Bessell’s filters 
llBessell. Castelli &: Plej|l998h are compared to those derived in this paper 
from Swift-WOT and XMM-Newton-OM observations of the PG 1553+113 
field. 


Band 

Star 1 

Star 2 

Star 3 

Star 4 

B 

14.503 (0.047) 

WEBT 
14.543 (0.050) 

16.399 (0.037) 

16.555 (0.039) 

V 

13.832 (0.027) 

13.923 (0.022) 

15.688 (0.017) 

15.771 (0.019) 

R 

13.465 (0.032) 

13.582 (0.029) 

15.277 (0.026) 

15.317 (0.029) 

I 

13.080 (0.056) 

13.230 (0.055) 

14.896 (0.051) 

14.893 (0.056) 

B 

14.497 (0.019) 

Swift-VNOT 

14.547 (0.019) 16.420 (0.051) 

16.580 (0.079) 

V 

13.885 (0.033) 

13.954 (0.030) 

15.721 (0.034) 

15.837 (0.043) 

B 

14.556 (0.022) 

XMM-Newton-OM. 

14.605 (0.022) 16.501 (0.023) 

16.652 (0.023) 

V 

13.838 (0.022) 

13.904 (0.022) 

15.676 (0.023) 

15.773 (0.023) 


1 


• 3 


• I 


.4 


Figure 1. Finding chart for PG 1553+113 with the stars of the 
optical photometric sequence labelled from 1 to 4. North is up 
and East on the left. The field of view is ~ 7 arcmin X 7 arcmin. 


were obtained after correcting for evident outliers, which 
are specially recognizable through a comparison among the 
source behaviour in different bands. Moreover, binning was 
used to reduce the noise of data acquired close in time by 
the same telescope. After the cleaning procedure, we were 
left with 3336 data points. No shift was given to correct for 
offsets between different data sets, since in the few cases 
where an offset was clearly identified, it was always lower 
than 0.05 mag. The final light curves are shown in Fig. 
where different symbols and colours highlight data from the 
various telescopes. The epochs of MAGIC, XMM-Newton 
and Swift observations are also indicated. 

Starting from the first data in January, the source 


brightness smoothly decreased by AR ~ 0.9 mag until June, 
when it reached a minimum; then it started to grow slowly. 
In the other, less-sampled optical bands the overall variation 
is about 0.5 mag. No significant flickering (fast continuous 
variability of the order of few tenths of mag on daily time 
scales) is visible. The three MAGIC observations occurred 
during the dimming phase, the last one very close to the 
minimum. In contrast, the XMM-Newton pointing met the 
source in the brightening phase. Swift observations are more 
distributed in time. 

Near-IR observations were performed at the Campo Im- 
peratore and Teide observatories (see Table [T]). The data 
sampling is not as dense as the optical one, but we can rec¬ 
ognize a similar trend. 

Radio data were acquired at 37 GHz with the 14 m an¬ 
tenna of the Metsahovi Observatory and at 43 and 8 GHz 
with the 32 m telescopes of the Noto and Medicina stations 
of the Istituto Nazionale di Astrofisica (INAF), respectively. 
The source is quite faint at radio wavelengths, its flux oscil¬ 
lating around 0.2 Jy. 


2.1 Colour analysis 

Although in the considered period the source was not partic¬ 
ularly active, we investigated the presence of possible spec¬ 
tral variations. We built B — R colour indices by coupling 
data acquired from the same telescope within 15 min. To in¬ 
crease the accuracy of the results, we only considered R- and 
H-band data with errors less than 0.02 and 0.04 mag, respec¬ 
tively, obtaining 169 values. The average is < B — R >=0.72 
and the standard deviation 0.03. The top panel of Fig. [3] 
highlights the points in the R-band light curve that have 
been used to build the B — R colour indices. These are plot¬ 
ted as a function of time in the middle panel, which also 
shows a smoothed B — R trend compared to a smoothed R- 
mag trencQ. The comparison suggests that most of the time 


^ Smoothing was obtained with a boxcar average of width 8, see 
http: //www.exelisvis.com/docs/SMOOTH.html 
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Figure 2. Cleaned light curves of PG 1553+113 built with data from the WEBT collaboration. Different colours and symbols highlight 
data from different telescopes. The yellow and green stripes mark the times of MAGIC and XMM-Newton observations, respectively; 
dotted lines those of Swift pointings. In the bottom panel purple circles, green squares and blue diamonds refer to flux densities at 43, 
37 and 8 GHz, respectively. 


B — R values increase with decreasing flux. This is more evi¬ 
dent in the bottom panel, where the behaviour of the colour 
index with brightness is displayed: a weak (linear Pearson’s 
correlation coefficient of 0.53 and Spearman’s rank corre¬ 
lation coefficient of 0.52 with two-sided significance of its 
deviation from zero of 3.6 x 10“^®), but clear, bluer-when- 
brighter trend is recognizable, which is often found in B L 
Lacs (see e.g. Ilkeiiri et al.l 12011 : Wierzcholska et al. gOlSll . 
It is interesting to notice that Wierzcholska et al. 1 2015ll did 
not find a definite trend when analysing the colour index 
behaviour in 2007-2012, when the source was in a brighter 


state than in 2013. However, the authors noticed that both 
bluer-when-brighter and redder-when-brighter trends can be 
found on shorter time scales. Indeed, a redder-when-brighter 
tendency is visible in Fig. Rafter JD = 2456500. 

In Fig. 13 we also show the value of the optical spectral 
index = [(B — R) — {Ab — Ar) — 0.306J/0.413, assuming 
that Fu oc u~°‘ is the de-reddened flux densit}0. The Galac¬ 
tic extinction values in the B and R bands, Ab and Ar, are 


^ We adopted effective wavele ngth and absolute flux values from 
iBessell. Gastelli &: Pled dlQQSh . 
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Figure 3. Top: the R-band light curve (black plus signs) where 
data points used to build the B — R colour indices have been 
highlighted with blue dots. Middle: B — R versus time; the red line 
represents the smoothed trend of the colour index; the green line 
the smoothed R-band light curve (R/2 —6.4). Bottom: B — R and 
optical spectral index Oo versus brightness; the red line represents 
a linear fit to the data. Horizontal dashed lines mark mean values. 


reported in Tabled We obtain a minimum and a maximum 
Qo of 0.60 and 0.97, with a mean index of < cto >= 0.80 and 
a standard deviation of 0.08. The spectral trend visible in 
the bottom panel then indicates that the synchrotron bump 
in the SED tends to peak very close to the optical band in 
faint states {uF^ oc 17““+^ ~ const), while in brighter states 
the peak shifts to higher frequencies (—a + 1 > 0). 

Using the same procedure we estimated the near- 
infrared spectral index from the J — K colour index as 
anir = [(J - K) - {Aj - Ak) - 0.987]/0.635. The aver¬ 
age value of 45 indices is < Qnir >= 0.67 with a standard 
deviation of 0.02. In this case there is no recognizable trend 
with brightness. 

We finally calculate the radio-optical spectral index 
Qro by considering the Medicina flux densities at 8 GHz, 
Fs, and simultaneous R-band data, a^o = 0.210 logTs + 
0.084 {R — An) — 0.731. We obtain values in the range 0.30- 
0.3d^, in agreemen t with the typical indices of HBL (e.g. 
iDoimto et al.llioolh . 


® We neglect fc-correction because we are dealing with an ob¬ 
ject of uncertain redshift. However, a roug h estimate assumin g 
2 = 0.5 and adopting the spectral indices of iFossati et al.l lllOOSh . 
indicates that the fc-correction should increase the index by only 
about 0.02. 


Table 3. Galactic extinction values (mag) used in this paper 
for the different photomet r ic ba nds. They are based on the 
ICardelli. Clayton fc Mathid lll989li laws with Ry = Ay jEiB — 
U) = 3.1, the standard value for the diffuse interstella r medi um, 
and on Ag = 0.224 from lSchlegel, Finkbeiner fc Dav3 lll99Sl 'l. In 
the case of UVOT and OM, the laws were convolved with the 
filter effective areas and source spectrum. 


Band 

Bessel 

UVOT 

OM 

W2 

- 

0.457 

0.484 

M2 

- 

0.481 

0.470 

Wl 

- 

0.394 

0.326 

U 

0.262 

0.275 

0.275 

B 

0.224 

0.229 

0.230 

V 

0.170 

0.173 

0.173 

R 

0.142 

- 

- 

I 

0.101 

- 

- 

J 

0.049 

- 

- 

H 

0.031 

- 

- 

K 

0.019 

- 

- 


Table 4. Swift pointings at PG 1553+113: starting times and 
duration of the XRT and UVOT exposures. 


ObslD 

Start time 

XRT (s) 

UVOT (s) 

00031368050 

2013-04-08 01:07:59 

1793 

1746 

00031368051 

2013-04-11 02:57:59 

2196 

2222 

00031368052 

2013-05-04 08:09:59 

1170 

1144 

00031368053 

2013-05-08 02:01:59 

1954 

1906 

00031368054 

2013-05-10 05:37:59 

859 

831 

00031368055 

2013-05-11 02:19:59 

1307 

1284 

00031368056 

2013-05-16 06:58:58 

1029 

1003 

00031368057 

2013-06-03 06:06:59 

1024 

997 

00031368058 

2013-06-08 06:11:59 

950 

920 

00031368059 

2013-06-08 22:33:59 

1040 

1011 

00031368060 

2013-06-09 00:13:58 

580 

551 

00031368061 

2013-06-11 21:13:59 

1974 

1918 

00031368062 

2013-06-13 05:06:59 

1014 

988 

00031368063 

2013-07-01 05:24:59 

1049 

1020 

00031368064 

2013-07-04 05:26:59 

1030 

1001 

00031368065 

2013-07-07 04:07:58 

1080 

1053 


3 OBSERVATIONS BY SWIFT 

During the WEBT campaign, the source was observed by 
the Swift satellite during 16 epochs. These are indicated in 
Fig. (2] Details on the observations are reported in Table |4] 
We rednced the data with the HEASOFT package ver¬ 
sion 6.15.1. 


3.1 UVOT 

The Swift sate llite carries a 30-cm Ultraviolet/Optical tele¬ 
scope (UVOT: lRoming et al.ll2005l ) equipped with v, b and u 
optical filters, and wl, m2 and w2 UV filters. To process the 
data, we used the UVOT calibration release 20130118 of the 
CALDB data base available at NASA’s High Energy Astro¬ 
physics Science Archive Research Center (HEASARC)Q. For 
each epoch, possible mnltiple images in the same filter were 

^ http://heasarc.nasa.gov/ 
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first summed with the task uvotimsum and then aperture 
photometry was performed with uvotsource. We extracted 
source counts from a circular region with 5 arcsec radius 
centred on the source and background counts from a circle 
with 20 arcsec radius in a source-free field region. 

The UVOT light curves are shown in Fig. H] On the 
left we plotted the observed magnitudes, while on the right 
we show flux densities after recalibration and correction for 
the Galactic extinction. Recalibration is necessary when a 
source has a spectral shape that differs from that of the ob¬ 
jects used to calibrate the instrument. The average b — v oi 
PG 1553-1-133 is 0.33, which is out of the colour range inside 
which the avail able UV count rate to flux conversion fa ctors 
(cfs) are valid dPoole et al.ll2008l : iBreeveld et al.ll2011^ . We 
thus interpolated the source spectrum with a log-linear fit 
and convolved it wi th the filter effecti ve areas to calculate 
new factors (see e.g. lRaiteri et al1l2010l i. The procedure was 
iterated to check the stability of the results. The ra t io be - 
tween these new cfs and the old cfs bv IBreeveld et aP (l201ll i 
is 0.997, 1.001, 1.015, 1.064, 0.995 and 1.014 in the v, b, u, 
wl, m2 and w2 bands, respectively. The most important ef¬ 
fect of the recalibration is thus to enhance the source flux 
density in the wl band by about 6%. 

Galactic extinction in eac h UVOT band was similarly 
estima ted by convolving the ICardelli. Glavton fc Mathid 
lll989ll laws with the filter effective areas and source spec¬ 
trum, adopting Rv = AyjEiJi — R) = 3.1 (the stan¬ 
dard value for the diffuse interstellar medium), and start¬ 
ing fr om an extinction of 0.224 mag in the Johnson’s B 
band (ISchlegel. Finkbeiner fc Davisl[l9^ '). The results after 
iteration are reported in Table (3] A discussion on the un¬ 
certainties affecting estimates of Gala ctic extinction, which 
are la rger in the UV, can be found in iFitzpatrick fc Massal 
|20 o 3). In the case of PG 1553-1-113, which has a low 
E{B — V) — 0.054 (see Table [3|, we would expect an er¬ 
ror of 0.02-0.04 mag on the UV de-reddened continuum due 
to uncertainties in the extinction alone. 

The UVOT light curves plotted in Fig. U confirm the 
trend shown by the ground-based optical (and near-IR) 
curves in Fig. [2l They also suggest an increasing variabil¬ 
ity amplitud e with frequency, as usually observed in BL Lac 
objects fe.g. iRaiteri et al .11201 Ol l. 

We also derived UVOT v and b magnitudes for the refer¬ 
ence stars and transformed them into Jo hnson’s magnitude s 
through the colour corrections given bv IPoole et al.l (l2008lj . 
The resulting V and B magnitudes, which differ from the 
V and b values by 0.009-0.017 mag, are reported in Table 
[21 They are in excellent agreement with the magnitudes de¬ 
rived from the WEBT data in the B band, while the UVOT 
data are a bit fainter in the V band. 


3.2 XRT 

The Swift X-reLy Telescope fXRT: [Burrows et al.ll2005h is a 
CCD imaging spectrometer. We reduced the data with the 
CALDB calibration files updated 20140120. We run the task 
xrtpipeline with standard screening criteria on the obser¬ 
vations performed in pointing mode. We found four obser¬ 
vations in Photon Counting (PC) mode and 16 observations 
in Windowed Timing (WT) mode, three of which were dis¬ 
carded because they had exposure times less than 30 s. The 
count rate ranges from 1.3 to 3.0 counts s“^ hence the four 
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Figure 4. UV—optical light curves of PG 1553-1-113 obtained from 
Swift-UYOTi data. Observed magnitudes are plotted on the left; 
maximum and minimum values are indicated in the lower left. 
The flux densities on the right were obtained after recalibration 
and correction for the Galactic extinction as explained in the text; 
maximum and minimum values are printed in the upper right. 


PC exposures need correction for pile-up. We analysed the 
PC images with the task ximage and modelled the wings of 
the source point spread function (PSF) with t he King’s func¬ 
tion representing the expected PSF of XRT l|Moretti et all 
I2OO5II . The extrapolation of the fit to the inner region al¬ 
lowed us to define the radius within which pile-up is impor¬ 
tant. This turned out to be 10 arcsec. We thus extracted 
source counts from an annular region with inner and outer 
radius of 10 and 75 arcsec, respectively, while we took the 
background from an annulus between 100 and 150 arcsec. As 
for the observations performed in WT mode, we extracted 
the source counts in a circular region with 70 arcsec radius, 
and the background counts from a region of the same size 
shifted along the window, away from the source PSF. The 
task xrtmkarf was used to create ancillary response files to 
correct for hot columns, bad pixels, and the loss of counts 
caused by using an annular extraction region in the pile-up 
case. 

We grouped each spectrum with the corresponding 
background, redistribution matrix (rmf), and ancillary (arf) 
files with the task grppha, setting a binning of at least 25 
counts for each spectral channel in order to use the 
statistic. The spectra were analysed with XSPEC version 
12.8.1. We adopted a Galactic ab sorption of Vh = 0.3 72 x 
10^^ cm“^ from the LAB sur vey dKalberla et al.l[2005l ) and 
the IWilms. Allen fc McCra\1 (l20od l elemental abundances. 
The spectra were fitted with both an absorbed power law 
N{E) = No E~^, where No represents the number of pho- 
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tons keV ^ cm ^ s ^ at 1 keV, and absorbed curved models. 
Application of a broken power law model 


N{E) = 


No 


E < Eo 
E ^ Eh 


led to large errors on the parameters or t o multiple solutions . 
We also tried the lo g-parabola model (iLandau et al.|[l986l : 
iMassaro et al.l[^004l ') 

N{E) = No 


(where Es is a scale parameter that we fixed equal to 1 keV) 
which has largely been used to fit the X-ray spectrum of this 
source (see Sect. li). 

Details of the power-law fits are given in Tabled where 
Column 1 gives the observation ID, Column 2 the XRT ob¬ 
serving mode, Column 3 the half-exposure JD, Column 4 
the 1 keV flux density corrected for the Galactic extinction. 
Column 5 the power-law photon index F, Column 6 the re¬ 
duced x^i and Column 7 the number of degrees of freedom. 
The average photon index is F = 2.14, with a standard de¬ 
viation of 0.06, and there is no significant correlation with 
brightness. Indeed, the linear Pearson correlation coefficient 
is only 0.46 and the variance is less than the mean square 
uncertainty, which means that possible variability is hidden 
by errors. Table [S] gives the same information as Table [S] for 
the log-parabolic fits, where the photon index F is replaced 
by the two parameters a and p. In the last Column we also 
report the T-test probability Pp (see below). 

Most of the power-law fits are statistically good (x?ed ~ 
1), with a few exceptions. The most critical case is the ob¬ 
servation at JD = 2456422.74, for which the fit cannot be 
improved by the curved model. In contrast, the data ac¬ 
quired on JD = 2456428.80 seem to be better reproduced 
by a log-parabola. To verify when the log-parabola offers a 
statistically better fit to the data than the power-law model, 
we applied the F statistic, which c an be used to test for 
an additional term in a model ('e.g.. I^otassov et al.l[200^ : 
lOrlandini et al.ll2012h . A low T-test probability Pp (e.g. be¬ 
low ~ 0.01) means that the simpler model is unlikely to be 
correct. The Pp in Table [6] covers a wide range of values, 
favouring the log-parabola at some epochs and disfavouring 
it in the others. The curvature parameter /3 corresponding 
to the cases with Pp < 0.01 goes from 0.27 to 0.53, indi¬ 
cating noticeable curvature, as expected. When instead /3 
approaches zero, Pp approaches one, or is not even com¬ 
putable. However, the large uncertainties on as well as its 
strong erratic variations prevent a reliable characterization 
of the X-ray spectral curvature with XRT. 

In Fig. [5] the de-absorbed X-ray light curve (obtained 
through power law fits) is compared with those at lower 
frequencies, showing the same decreasing trend of the near- 
IR to UV light curves. 


4 OBSERVATIONS BY XMM-NEWTON 

A long calibration pointing at PG 1553-1-113 was per¬ 
formed by the XMM-Newton satellite on 2013 July 24-25 
(JD = 2456498.13560-2456498.52237). The total exposure 


2013.25 2013.30 2013.35 2013.40 2013.45 



6380 6400 6420 6440 6460 

Julian Date - 2450000 


Figure 5. Multiwavelength light curves of PG 15534-113 around 
the period of the MAGIC pointings (yellow stripes). In the top 
panel the de-absorbed flux density at 1 keV is given, separating 
the data obtained in PC (red dots) from WT (blue dots) mode. 
Flux densities in the uil, J? and H bands are corrected for Galactic 
extinction as explained in the text. Different colours and symbols 
in the bottom panel, showing the radio flux densities, have the 
same meaning as in Fig. [2] 


time was 34500 s. We processed the data with the Science 
Analysis Softwar^fl (SAS) package version 13.5.0. 


4.1 OM 

The Optical Monitor llMason et al.l 1200 ll ) onboard XMM- 
Newton performed nine exposures, six of which with the 
optical and UV wide-band filters, two with the UV grism 
(grisml), and one with the optical grism (grism2). We ran 
the pipeline omichain and analysed the V, B, Wl, M2 and 
W2 images with the omsource task. The resulting source 
magnitudes, already adjusted to the Johnson’s system, are 
reported in Table [T] The uncertainties take into account 
both the measurement and the calibration errors. The mea¬ 
surement error is very small, of the order of 0.001-0.003 mag, 
in all bands but W2, where it is 0.018 mag. The errors in the 
calibration are estimated to be 2% in V and B bands and 
10 % in U bantfl. We assumed a 10% error also in the UV 
bands. We notice that the OM V and B magnitudes fairly 

® http://xmm.esac.esa.int/sas/current/documentation/ 

® http://xmm2.esac.esa.int/docs/documents/CAL-TN-0019.pdf 
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Table 5. Results of fitting the Swift-X.KT spectra with an absorbed power-law model with photon index F and fixed Galactic absorption 
A^h = 0.372 X 10^1 cm-2. 


ID 

Mode 

JD 

.F’lkeV [I'Jy] 


r 

X?ed 

d.o.f. 

00031368050 

PC 

2456390.55759 

2013 

11.46 (11.01-11.91) 

2.14 

(2.09-2.20) 

0.95 

77 

00031368051 

PC 

2456393.63631 

9.83 (9.40-10.25) 

2.12 

(2.05-2.18) 

1.17 

71 

00031368051 

WT 

2456393.63631 

11.30 (10.56-12.05) 

2.19 

(2.09-2.29) 

1.07 

46 

00031368052 

WT 

2456416.84703 

8.01 (7.67-8.34) 

2.15 

(2.09-2.22) 

0.99 

92 

00031368053 

PC 

2456420.59602 

8.89 (8.51-9.26) 

2.15 

(2.09-2.20) 

1.07 

72 

00031368054 

WT 

2456422.73968 

9.48 (9.06-9.89) 

2.14 

(2.08-2.21) 

1.40 

81 

00031368055 

PC 

2456423.60478 

7.89 (7.46-8.32) 

2.11 

(2.03-2.19) 

0.97 

45 

00031368056 

WT 

2456428.79691 

8.74 (8.36-9.11) 

2.19 

(2.13-2.25) 

1.21 

86 

00031368057 

WT 

2456446.76078 

7.56 (7.20-7.93) 

2.09 

(2.02-2.16) 

1.09 

72 

00031368058 

WT 

2456451.76382 

6.58 (6.23-6.94) 

2.08 

(2.00-2.16) 

1.10 

63 

00031368059 

WT 

2456452.44628 

7.29 (6.95-7.62) 

2.24 

(2.17-2.32) 

0.85 

74 

00031368060 

WT 

2456452.51305 

6.74 (6.28-7.20) 

2.01 

(1.92-2.11) 

0.93 

41 

00031368061 

WT 

2456455.39613 

7.01 (6.77-7.26) 

2.17 

(2.12-2.22) 

0.98 

131 

00031368062 

WT 

2456456.71905 

6.70 (6.37-7.03) 

2.11 

(2.04-2.19) 

0.93 

73 

00031368063 

WT 

2456474.73176 

9.38 (8.99-9.78) 

2.15 

(2.09-2.21) 

1.17 

91 

00031368064 

WT 

2456477.73303 

10.82 (10.37-11.26) 

2.20 

(2.14-2.26) 

0.92 

94 

00031368065 

WT 

2456480.67845 

10.35 (9.95-10.74) 

2.21 

(2.16-2.27) 

1.01 

102 

00031368009 

PC 

2455083.53235 

2009 

2.63 (2.54-2.72) 

2.37 (2.33-2.42) 

1.11 

92 

00031368010 

PC 

2455108.09667 

2.18 (2.08-2.28) 

2.43 (2.36-2.51) 

1.06 

54 

00031368035 

WT 

2456046.57654 

2012 

17.12 (16.61-17.63) 

2.16 (2.12-2.20) 

1.15 

144 


Table 6. Results of fitting the XRT spectra with an absorbed log-parabola model N{E) = No (E/Es) “ /31og(B/Bs) with Es = IkeV 
and fixed Galactic absorption Wh = 0.372 X 10^^ cm“^. 


ID 

Mode 

JD 

.P’lkeV [Z^Jy] 

a 

0 

Xred 

d.o.f. 


Pf 

00031368050 

PC 

2456390.55759 

2013 

12.26 (11.66-12.86) 

2.06 (1.99-2.13) 

0.37 (0.20-0.56) 

0.80 

76 

2 

X 10-* 

00031368051 

PC 

2456393.63631 

10.19 (9.66-10.72) 

2.06 (1.98-2.14) 

0.22 (0.04-0.42) 

1.14 

70 


0.09 

00031368051 

WT 

2456393.63631 

11.44 (10.54-12.31) 

2.16 (2.02-2.29) 

0.09 (-0.18-0.40) 

1.09 

45 


0.69 

00031368052 

WT 

2456416.84703 

8.17 (7.76-8.58) 

2.12 (2.03-2.20) 

0.14 (-0.06-0.35) 

0.98 

91 


0.17 

00031368053 

PC 

2456420.59602 

9.77 (9.26-10.27) 

2.03 (1.94-2.10) 

0.52 (0.32-0.73) 

0.79 

71 

2 

X 10-® 

00031368054 

WT 

2456422.73968 

9.56 (9.07-10.07) 

2.13 (2.04-2.21) 

0.06 (-0.13-0.26) 

1.42 

80 


- 

00031368055 

PC 

2456423.60478 

8.26 (7.71-8.81) 

2.05 (1.95-2.15) 

0.28 (0.03-0.54) 

0.92 

44 


0.07 

00031368056 

WT 

2456428.79691 

9.41 (8.94-9.88) 

2.06 (1.97-2.15) 

0.53 (0.32-0.75) 

1.01 

85 

6 

X 10“® 

00031368057 

WT 

2456446.76078 

7.96 (7.51-8.41) 

1.99 (1.88-2.08) 

0.37 (0.15-0.62) 

0.99 

71 

5 

X 10-® 

00031368058 

WT 

2456451.76382 

6.96 (6.53-7.36) 

1.97 (1.85-2.08) 

0.42 (0.15-0.72) 

1.01 

62 


0.01 

00031368059 

WT 

2456452.44628 

7.42 (6.96-7.83) 

2.22 (2.13-2.30) 

0.12 (-0.10-0.35) 

0.85 

73 


0.32 

00031368060 

WT 

2456452.51305 

7.23 (6.66-7.81) 

1.88 (1.74-2.02) 

0.53 (0.19-0.90) 

0.78 

40 

5 

X 10“® 

00031368061 

WT 

2456455.39613 

7.31 (7.00-7.61) 

2.10 (2.03-2.17) 

0.28 (0.12-0.45) 

0.93 

130 

5 

X 10-3 

00031368062 

WT 

2456456.71905 

6.83 (6.44-7.22) 

2.08 (1.97-2.17) 

0.14 (-0.08-0.37) 

0.92 

72 


0.18 

00031368063 

WT 

2456474.73176 

9.53 (9.05-10.01) 

2.12 (2.04-2.20) 

0.11 bo.09-0.31) 

1.18 

90 


0.63 

00031368064 

WT 

2456477.73303 

10.82 (10.29-11.35) 

2.20 (2.12-2.27) 

0.00 bo.18-0.19) 

0.93 

93 


- 

00031368065 

WT 

2456480.67845 

10.70 (10.22-11.19) 

2.15 (2.08-2.23) 

0.23 (0.06-0.42) 

0.97 

101 


0.02 

00031368009 

PC 

2455083.53235 

2009 

2.82 (2.70-2.95) 

2.32 (2.26-2.38) 

0.39 (0.22-0.57) 

0.95 

91 


10-"* 

00031368010 

PC 

2455108.09667 

2.25 (2.12-2.39) 

2.42 (2.34-2.50) 

0.17 (-0.04-0.41) 

1.04 

53 


0.16 

00031368035 

WT 

2456046.57654 

2012 

17.91 (17.27-18.56) 

2.10 (2.04-2.15) 

0.27 (0.14-0.41) 

1.07 

143 
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Table 7. Results of the XMM-Newton-OM photometry on PG 
1553+113. 


Filter 

Exp [s] 

Mag 



V 

1700 

14.561 (0.022) 

5.72 

6.60 

B 

2001 

14.986 (0.022) 

6.65 

5.55 

U 

2000 

14.09 (0.11) 

9.05 

4.60 

Wl 

3999 

13.84 (0.11) 

10.6 

4.04 

M2 

3999 

13.79 (0.11) 

13.6 

3.74 

W2 

4000 

13.87 (0.11) 

14.4 

3.36 


^ Observed flux densities in units 10 erg cm ^ s ^ A ^ 
^ De-absorbed flux densities in mjy 


match the simultaneous values by the WEBT shown in Fig. 

m 

We also extracted the V and B magnitudes of the ref¬ 
erence stars, to compare their space photometry with the 
ground-based one. As can be seen in Tabled the OM values 
in V band are in excellent agreement with those obtained 
with the WEBT data, with deviations of less than 0.02 mag. 
The difference in B band is more pronounced, up to ~ 0.1 
mag. 

The OM grisml and grism2 provide low-resolution spec¬ 
tra in the range 1800 to 3600 A and 3000 to 6000 A, respec¬ 
tively. Exposure times were 4400 and 4260 s for the two 
grisml observations and 4000 s for the grism2. We reduced 
these data with the task omgchain and checked the results 
with omgsource. The resulting UV spectra are very noisy 
and no reliable feature can be recognized. The optical spec¬ 
trum is contaminated by zero-order features. 

The grism spectra are shown in Fig. [6l where they are 
compared with the OM broad-band photometry. Flux den¬ 
sities in the V, B, U, Wl, M2 and W2 bands were obtained 
by multiplying the count rates for the conversion factors 
derived from observations of white dwarf standard starj^. 
They are also reported in Table [T] We assumed a conser¬ 
vative 10% error on these values. There is a fair agreement 
between the results of the OM spectroscopy and photometry. 

In Table [7] we also give de-absorbed flux densities in 
mJy. The values of Galactic extinction for the OM fil- 
ters calculated by integrating the source spectrum with the 
ICardclli. Clayton fc Mathis! ( 19891 ') laws and filter effective 
areas are reported in Table [3] 

4.2 EPIC 

The European Photon Imaging Camera (EPIC) on¬ 
board XMM-Newton carries three detectors: MOSl, MOS2 
llTurner et al.l l20()i] ~l and pn (IStrtider et al.lliooH l. Because 
of the expected high source brightness, exposures were per¬ 
formed in small window mode with medium filter. We pro¬ 
cessed the data with the emproc and epproc tasks of the SAS 
package. Possible high-background periods were removed by 
asking that the count rate of high-energy events (> lOkeV) 
was less than 0.35 and 0.40 cts s“^ on the MOS and pn de¬ 
tectors, respectively. Out-of-time events are visible in the pn 

http://xmm.esac.esa.int/sas/8.0.0/watchout/. We notice that 
when using conversion factors based on the Vega spectrum, a dip 
in U band appears. 



Figure 6. XMM-Newton-OXl grism spectra and broad-band pho¬ 
tometry. The green and red lines represent the grisml (UV) spec¬ 
tra, while the blue line shows the grism2 (optical) spectrum. Black 
dots mark the flux densities in the V, B, U, Wl, M2 and W2 
Alters. A conservative 10% error is assumed on the flux. 


image, but because of the small window mode, their frac¬ 
tion is ~ 1% only and a correction is not necessary. We 
extracted source counts from a circular region with 40 arc- 
sec radius. In the case of MOS detectors, background was 
extracted from a circle of 80 arcsec radius on an external 
CCD, while in the case of pn it was taken from a circle of 
40 arcsec radius on the same CCD, as far as possible from 
the source. We selected the best calibrated single and dou¬ 
ble events only (PATTERN<=4), and rejected events next 
to either the edges of the CCDs or bad pixels (FLAG==0). 
The absence of pile-up was verified with the epatplot task. 

As in the XRT data analysis, we grouped each spectrum 
with the corresponding background, redistribution matrix 
(rmf), and ancillary (arf) files with the task grppha, set¬ 
ting a binning of at least 25 counts for each spectral chan¬ 
nel in order to use the statistic. The three spectra were 
analysed with XSPEC version 12.8.1. We adopted a Galactic 
abso rption of Ah = 0.372 x 10^^ cm“^ from the LAB sur- 
yev llKalberla et al.ll200^ and the lWilms. Allen fc McCravl 
ll200oh elemental abundances. We fitted the three EPIC 
spectra together. An absorbed power law gives a poor fit to 
the spectrum, even when letting the absorption free. The X- 
ray spectra are better reproduced by absorbed curved mod¬ 
els such as a broken power law or a log-parabola (see the 
previous section). Table [S] shows the results of the spec¬ 
tral fitting with the different models. Both curved mod¬ 
els give a reasonably good fit according to the Xred value. 
Unfortunately, this is not a case wh ere we can use the F- 
statistic to see which model is better llProtassov et al1l2002l : 
lOrlandini et al.ll2012ll . so we consider both as reliable. 

In the parabolic case, the observed fluxes in the 0.3- 
10, 0.3-2 and 2-10 keV energy ranges are (6.64 ± 0.02), 
(4.27 ± 0.01) and (2.37± 0.01) x 10~^^ ergcm“^ s“^, respec¬ 
tively. The corresponding de-absorbed values are 7.62, 5.25 
and 2.37 x 10“^^ ergcrn”^ s“^. The three EPIC spectra with 
folded model and data-to-model ratio are shown in Fig. [7] 

We finally extracted source X-ray light curves with dif¬ 
ferent time binning, but they did not show significant vari¬ 
ability. 
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Table 8. Results of fitting the XMM-Newton-'EiPlQ spectra with different models: power law (pow), broken power law (bknpow) 
and log-parabola (logpar). For the broken power law and log-parabola models the absorption was fixed to the Galactic value A^h = 
0.372 X 10^^ cm“^, while for the single power law models both free and fixed absorption were tried. 


Model 

pow 

ElkeV [hJy] 

10.48 ± 0.02 

Ah [10^1 cm-2] 
0.372 

r 

2.319 ±0.002 


X?ed 

1.56 

d.o.f. 

3270 

Model 

pow 

FlkeV [^Jy] 

11.15 ±0.04 

Ah [10^1 cm-2] 
0.550 ±0.009 

r 

2.395 ± 0.004 


Ared 

1.19 

d.o.f. 

3269 

Model 

bknpow 

FlkeV [^Jy] 

10.75 ±0.03 

Ti 

2.246 ± 0.005 

Vi 

2.464 ±0.010 

E6 [keV] 
1.62 ±0.06 

^red 

1.08 

d.o.f. 

3268 

Model 

logpar 

FlkeV [^Jy] 

10.83 ±0.03 

OL 

2.283 zb 0.003 

d 

0.166 ±0.007 

Es [keV] 

1 

Ared 

1.09 

d.o.f. 

3269 



. . ' ... ,"||M ; , 

0.5 1 2 5 10 

Energy {keV) 


Figure 7. X-ray spectra of PG 1553-1-113 acquired by the EPIC 
instrument onboard XMM-Newton on 2013 July 24—25. Black, 
red and green symbols represent MOSl, MOS2 and pn data, 
respectively. The folded model, an absorbed log-parabola with 

= 0.372 X 10^^cm“^, is shown by solid lines of the same 
colour. The ratio between the data and the folded model is plot¬ 
ted in the bottom panel. 

5 SPECTRAL ENERGY DISTRIBUTION 

In the following, we present SEDs of the source from the 
near-infrared to the X-ray energies, i.e., the frequency range 
where the synchrotron bump achieves its maximum in the 
high-energy peaked BL Lac objects such as PG 1553-1-113. 
For the UV data of both UVOT and OM we assumed a 10% 
error on the flux. The errors on the WEBT data also include 
uncertainties on the calibration. The S'tui/t-XRT spectra are 
here modelled with power laws. 

At the epoch of the first MAG/C pointing, from JD = 
2456390.6326 to JD = 2456393.7306, two Swift observations 
were performed: the Erst one just before the start and the 
second by the end of the MAGIC pointing. Figure |8] shows 
the corresponding SEDs of the source; X-ray and UV data 
from Swift are complemented by simultaneous optical and 
near-IR data acquired by the WEBT. In particular, in the 
hrst SED (red) the J, H and K data are from Campo Im- 
peratore, while the R-band datum is the average of eight 
cleaned data from Abastumani. Both a PC-mode and a WT- 
mode spectra were available for the second SED (blue); we 
chose to plot the results of the PC mode because of the larger 
number of counts. Near-IR points for this epoch represent 


averages of simultaneous WEBT data from Campo Impera- 
tore and Teide; optical B, V, R and I data are mean values 
from simultaneous observations at Rozhen and a dense mon¬ 
itoring at Mt. Maidanak. We notice that the WEBT points 
nicely continue the spectral shape traced by the optical- 
UV data from Swift-IJYOT and that the U-band points 
from UVOT and from the WEBT available in the second 
epoch are in excellent agreement. The only critical point is 
in the B band, where the WEBT fluxes are a bit fainter than 
the UVO T values. This prob lem has already been found for 
3C 454.3 llRaiteri et al.ll201l| j. 

In the period of the second MAGIC pointing (JD = 
2456420.5826-2456423.6299) there were three Sm/t observa¬ 
tions. WEBT BVRI data corresponding to the first epoch 
are from the Crimean Observatory, complemented by ob¬ 
servations at the Abastumani Observatory in the R band. 
WEBT data for the second pointing were provided by the 
Crimean, Mt. Maidanak and New Mexico Skies observato¬ 
ries, while for the third pointing WEBT data were acquired 
at the Valle d’Aosta Observatory. 

The third MAG/Gpointing, from JD = 2456452.4486 to 
JD = 2456455.5083, saw again three Swift pointings. In the 
figure, we only plot the first and last ones, since the second 
is very close in time to the first but lacks observations in 
the m2 and v bands. WEBT data for these two epochs are 
from the Mt. Maidanak Observatory, complemented in the 
last epoch by observations from the Crimean Observatory. 
No near-IR data are available during the second and third 
MAGIC observations. 

A log-parabo lic model llLandau et al.l Il986l : 
iMassaro et al.l |2004^ applied to the near-IR-X-ray SEDs 
gave a poor fit to the X-ray data (x?ed = 7-14). Hence 
we htted all SEDs with a log-cubic model (x?ed — 2-7), 
which suggests synchrotron peak frequencies in the range 
log i/p = 15.3-15.6, i.e. Vp = 8-16 eV, and a general increase 
with the source brightness. An increase of the synchrotron 
peak frequency with flux has already been found for various 
other blazars, notabl y for another hig h-energy peaked BL 
Lac object: Mkn 501 (IPian et al.lll99!J l. Note however that 
the fit cannot account for the downturn in the UV band. 

The SED corresponding to the XMM-Newton pointing 
of 2013 July 24-25 (JD = 2456498.13560-2456498.52237) is 
shown in Fig. Both the log-parabola (blue) and the bro¬ 
ken power-law (green) hts to the EPIC spectra are plotted. 
The OM spectrum is hard. The simultaneous WEBT data 
points are from the Crimean Observatory. The B and V 
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ground-based fluxes appear in agreement, within the errors, 
with those derived from the space observations. The dotted 
and dashed lines represent log-parabolic and log-cubic fits to 
the data, respectively^. They highlight that the OM UVM2 
and UVVT2 data points are a bit too high to smoothly con¬ 
nect with the X-ray spectrum, whose curvature cannot be 
fully accounted for. The fits suggest a synchrotron peak lo¬ 
cated at log r'p = 15.8-15.9, i.e. Up ~ 26-33 eV. For com¬ 
parison, a parabolic fit through the SED including XMM- 
Newton observations in 2001, when the source was 1.5 times 
brighter, gave logr'p = 15.7 dPerlman et al.ir2005h . In con¬ 
trast, the log-parabolic curvature of the 2013 X-ray spec¬ 
trum alone would imply logj/p ~ 16.5 at logvFi, —10.46, 
incompatible with a single synchrotron peak that also in¬ 
cludes the UV and optical data. 


6 DISCUSSION 

In the previous section we presented broad-band SEDs of 
PG 1553-1-113 at different epochs, where UV and X-ray data 
from either Swift or XMM-Newton were available. The SEDs 
built with Swift data present a convexity in the optical-UV 
that makes a smooth connection of the UV with the X-ray 
spectra hard to trace, especially when adopting a curved 
model for the X-ray spectrum. This problem does not de¬ 
pend on the recalibration procedure we applied, since stan¬ 
dard calibration even worsens the picture. The errors on the 
UVOT points however are rather large as well as the uncer¬ 
tainty affecting the XRT spectra. 

In contrast, the SED obtained with the XMM-Newton 
data is very precise and the problem of connecting the UV to 
the X-ray spectrum is mitigated, but not completely solved 
because of the curvature inferred from the spectral analysis 
of the EPIC data. The polynomial fits suggest a synchrotron 
peak at a higher frequency than in the Swift case. It is in¬ 
teresting to notice that a simil ar X-ray spectral shape was 
found by IPerlman et al.l (( 200 I) while analysing an XMM- 
Newton observation of September 2001 (a = 2.331 q qj, P = 
0-13tn'n^), and an even more p ronounced curvature was ob¬ 
tained by iReimer et al.l (l2008l l from a Suzaku observation 
in July 2006 (a = 2.19 ± 0.01, /3 = 0.26 ± 0.01). In both 
epochs the X-ray flux was about 1.5 times higher than in 
the XMM-Newton observation of 2013. Noti c eable curva¬ 
tures were also estimated by iMassaro et al.l (l2008fl from 
Swift observations in 2005 (a = 2.11-2.21, /? = 0.23- 
0.36) and especially from a BeppoSAX pointing in 1998 
(a = 2.17±0.07, P = 0.63±0.08). In all these cases, the log- 
parabolic fit of the X-ray spectrum gave a synchrotron peak 
energy of about 0.5-0.7 keV, i.e. more than 10 times higher 
than we found by fitting the SED over a broader energy 
range. Physical explanations of the X-ray spectral curva¬ 
ture o bserved in BL Lac objects were given bv iMassaro et all 
ll2004h by assuming that the probability for a particle to in¬ 
crease i ts energy is a d e creasi ng function of the energy itself, 
and by I Perlman et aP ([200^ as due to episodic particle ac¬ 
celeration. In the latter model particle acceleration occurs 
with a typical timescale T and the maximum synchrotron 


The E-test probability is 0.17, which does not clearly favour 
any of the two models. 
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Figure 8. Spectral energy distribution of PG 1553+113 during 
the first (top), second (middle) and third (bottom) MAG/C obser¬ 
vations. Data from the WEBT are displayed as diamonds, while 
those from Swift-WOT as circles. The errors on the WEBT data 
take into account also the uncertainties on the calibration. The 
dashed lines represent log-cubic fits, which suggest a general shift 
of the synchrotron peak (marked by vertical dotted lines) towards 
higher frequencies when the source brightness increases. 


energy Ei is reached at time T, when the acceleration stops. 
The observed photon spectrum is given by: 




dE 


Ei/2 [eUL - E[''^) 


1 / 2 , 


where Umax is the maximum energy that would be reached 
under steady injection conditions. We show in Fig. [9] a fit of 
this model to the 2013 EPIC spectrum, normalized to 0.3 
keV and with the choice of parameters s = 3.0, Ei = 25 keV 
and Umax = 100 keV. The model gives a fair fit of the EPIC 
spectrum, though predicts a slightly stronger curvature at 
high energies, while the data show a more pronounced bend¬ 
ing in soft X-rays. 

With the aim of clarifying the shape of the source spec- 
tra around the synchro tron peak we looked for other data. 
iFalomo fc Trevea (Il990fl analysed multiwavelength observa- 
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Figure 9. Spectral energy distribution of PG 1553+113 during 
the XMM-Newton pointing of 2013 July 24-25 (JD = 2456498). 
Data from the WEST are displayed as diamonds, while those 
from XMM-Newton-OM as squares. The errors on the WEBT 
data take into account also the uncertainties on the calibration. 
The EPIC spectra include errors on the free parameters of the log- 
parabola (blue) and broken power law (green) models. The dotted 
and dashed lines are log-parabolic and log-cubic fits, respectively, 
suggesting a synchrotron peak at about 26-33 eV. The red dot- 
dashed line represent s the episodic particle acceleration model by 
IPerlman et all ll2n05h . 

tions of PG 1553+113 including four pointings by the lUE 
in 1982-1988. From their paper we extracted the infrared- 
to-UV SED corresponding to a high state of the source in 
August 1988. As can be seen in Fig. 1101 it shows a curvature, 
suggesting a synchrotron peak around log t'p = 14.9. 

In the same figure we plotted the HST-COS 
spectrum acquired on 2009 September 22 and down¬ 
loaded from the Space Telescope Science Institute 
arc hivJ^. The spectrum was de-ab sorbed according to 
the lCardelli. Clayton fc Mathid (1 19891 ) law and smoothecF^. 
The wealth of int erstellar and intergalac tic medium features 
were analysed bv iDanforth et all (l20ld ). who used them to 
set constraints on the source redshift (see the Introduction). 
The HST spectrum covers the wavelength range 1120- 
1800 A and is soft. It refers to a fainter state of the source 
than the lUE spectrum. Swift pointed at the source on 2009 
September 9 and October 3, i.e. 13 days before and 11 days 
after the HST observations. We analysed the correspond¬ 
ing data following the same procedure described in Sect. [I] 
The results are shown in Fig. 1101 while details on the fit 
of the XRT spectra are given in Table [S] and [B] While the 
X-ray flux decreased by ~ 20% from September 9 to Octo¬ 
ber 3, the optical-UV fluxes are similar in the two epochs. 
Hence, if we assume that in that period the UV flux was 
stable, the composite Swift plus HST SED can be consid¬ 
ered as indicative of the source spectral shape at that time. 
A fourth-order polynomial fit shows how the X-ray spectra 
can smoothly connect with the optical-UV ones, suggesting 

http://mast.stsci.edu 

We also cut the deep absorption line at A ~ 1210-1220 A for 
sake of clarity. 


Figure 10. Spectral energy distributions of PG 1553+113 in dif¬ 
ferent brightness states. The Swift data of 2009 September 9, 2009 
October 3, 2012 April 29 and 2013 May 10 are shown in red, blue, 
orange and green, respectively. The X-ray spectra are modelled 
with power laws; in the high and low states also log-parabola fits 
are shown. The de-absorbed HST-COS spectrum of 2009 Septem¬ 
ber 22 is plotted in black and prolongs the Swift-VYOT spec¬ 
tra of 2009 September 9 and October 3. The dashed lines rep¬ 
resent fourth-order polynomial fits to the near-IR to power-law 
X-ray spectra (including the HST-COS data in 2009); the dot- 
dashed line is a fit to the optical to log-parabolic X-ray spec¬ 
trum of 2012 April 29. The purple crosses and lines refer to near- 
IR-UV data acquired i n August 1988 close to HIE observations 
llFalomo fc Trev5ll990l ). 

that the synchrotron peak is at log Vp « 15.0, i.e. about 4 
eV. This would further confirm the trend of decreasing Vp 
with decreasing X-ray flux already seen for the 2013 Swift 
observations. We finally notice that in 2009 the X-ray flux 
was 3-6 times lower than in 2013 and the X-ray spectrum 
definitely softer, while the optical-UV level was similar (see 
Fig.®. 

Hence, both the lUE and the HST-COS spectra con¬ 
firm the soft Swift-CYOT spectra, so we can guess that 
the XMM-Newton-OM data may suffer calibration prob¬ 
lems. Yet, a recalibration of the OM data following the same 
prescriptions that we used for UVOT does not produce ap¬ 
preciable differences. 

The 2013 and 2009 SEDs that we have analysed show 
faint states of the source. In order to investigate the SED 
behaviour at brighter levels, we analysed Swift data taken 
on 2012 April 29, a t the peak of the 2012 outburst (see also 
lAleksic et al.ll2015h . The results of spectral fitting are re¬ 
ported in Tables [B] and [B] while the corresponding SED is 
displayed in Fig. HOI The cubic fit suggests a synchrotron 
peak at log Vp « 15.3, i.e. about 8.3 eV. Because of the 
higher statistic compared to the other S'uii/f-XRT observa¬ 
tions, we also show the result of a log-parabolic fit on the 
X-ray spectrum. This implies a more pronounced curvature 
than in the XMM-Newton case and suggests an inflection 
point around log v = IQ. 

The problem of connecting the optical-UV spectrum 
with the X-ray one presents no simple solution. One may 
wonder whether t he choice of Galactic extinction can p lay a 
role. We followed I Schlegel. Finkbeiner fc David lll998l) : had 
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we used the recalibration bv ISchlaflv fc Finkbeineil ll201lli . 
which implies a ~ 20% lower extinction values, we would 
have obtained a small shift of the optical-UV SEDs down¬ 
ward, by ~ 0.01 in the V band, up to 0.04 in the w2 
band. This would slightly improve the match between the 
XMM-Newton UV and X-ray spectra, but would worsen it 
in the Swift case. In general, it seems hard to find a solu¬ 
tion that suits both the Swift and XMM-Newton data just 
playing with the Galactic extinction uncertainties (see also 
Sect. [6^ . 

Incidentally, we notice that similar connection problems 
were found for ot her BL Lac objects, both LBL such a s 
BL Lacertae itself llRavasio et al.l 20031: Raiteri et aLl l2009l l 
and HBL such as Mkn 421 1 Massaro et al.l 12004 ) and H 
1722+119 (Ahnen et al. 2015, in preparation). Moreover, 
inflection points in the infrared portion of the SED have 
been revealed in 3 C 454.3, PKS 1 510—089 and 3C 279 
dWehrle et al.l[^012l : iNalewaiko et al.]l2012l : iHavashida et al.) 
I 2 OI 2 I I. All these evidences suggest that a full understanding 
of blazar emission requires a more complex picture than that 
is commonly assumed. 


6.1 Helical jet interpretation 

In the following we will analyse the source SEDs in an inho¬ 
mogeneous curved jet scenario, where synchrotron radiation 
of decreasing frequency is emitted from jet regions at in¬ 
creasing distance from the jet apex, which have different 
viewing angles and thus different Doppler factors. A helical 
j et mo del of this kind was developed by IVillata fc Raiteril 
lll999lf and was successfully applied to interpret the ob- 


served broad-band flux 

variability in other BL Lac objects 

fe.g. iRaiteri et al.l 

1999 

lOstorero. Villata & Raiteril |2004|; 

iRaiteri et al.l 2009(. 

2010 

). According to this model, the jet 


viewing angle varies along the helical path as 


cosff(z) = COS')/) cos C + sin ^ sin C cos((f — az) , (1) 

assuming that the helical jet axis lies along the a coordi¬ 
nate of a 3-D reference frame. The angle ‘ip is defined by 
the helix axis and the line of sight, ^ is the pitch angle, 
(f) is the azimuthal difference between the line of sight and 
the initial direction of the helical path, and a is the “cur¬ 
vature”, defining the azimuthal angle ip{z) = az along the 
helical path. Each slice of the jet can radiate, in the plasma 
rest frame, synchrotron photons from a minimum frequency 
I'min f o a maximum one r'max • The observed flux density has 
a power-law dependence on the frequency (with spectral in¬ 
dex oo) and a cubic dependence on the Doppler beaming 
factor 5 = [rb(l — /?cos0)]“^, where jS is the bulk veloc¬ 
ity of the emitting plasma in units of the speed of light, 
Tb = (1 — the corresponding bulk Lorentz factor, 

and 6 is the viewing angle of equation [T] The variation of 
the viewing angle along the helical path implies a change of 
the beaming factor. As a consequence, the flux at v peaks 
when the part of the jet mostly contributing to it has mini¬ 
mum 9. 

In Fig. [ 11 ] we show fits to three different brightness 
states of the source obtained with this model, where only 
the three parameters defining the jet orientation (a, tp and 
(p) have been changed from one fit to the other, leaving the 
intrinsic flux unchanged. The high state is obtained with 
a = 45°, Ip = 25.7° and ip = 17.5°, the intermediate state 


with a = 48°, ^p = 24° and (p = 18°, and the low state with 
a = 70°, Ip = 24.8° and cp = 26.5°. As for the other param¬ 
eters introduced above, we fixed ( = 30°, log z.'(n;j,(0) = 16.7 
and logi^max(O) = 19-3 at jet apex {z = 0), ao = 0.5, and 
Eb = 10. We selected the model parameters in particular to 
obtain i) an X-ray spectrum with the same curvature as that 
derived from the XMM-Newton data, but ii) a UV spectral 
shape more similar to that given by the Swift-WOT data, 
iii) an inflection point in the SED of the low X-ray state 
that can reproduce the narrow synchrotron peak marked by 
the Swift-WOT and HST-COS data, and iv) sufficient flux 
in the radio band to fit the radio observations at 43 GHz 
performed close to the XMM-Newton pointing of 2013 July 


24-25. 

While we refer to IVillata fc Raiteril (| 19991 . and to the 
other papers mentioned at the beginning of this section) for 
a full mathematical description of the helical jet model, we 
show in Fig.[T2]the trends of some physical and geometrical 
quantities with the distance z along the helix axis. In the top 
panel one can see that the maximum and minimum emitted 
frequencies decrease along the jet, while the emissivity first 
rapidly increases and then slowly decrease£3 In the bottom 
panel we compare the behaviour of the viewing angle 9 and 
Doppler factor S for the three fits to the source SEDs pre¬ 
sented in Fig. 111! Close to the jet apex, where the highest 
frequencies are emitted, the Doppler factor is greater for the 
high-state fit and lower for the low X-ray state, and this is 
also true, with different amplitude, along all the jet. For all 
cases the maximum beaming occurs around z = 0.38-0.4, 
the region from which frequencies in the range 10^° to 
lO^'* Hz are emitted. The peak of the Doppler factor is lower 
in the intermediate state than in the low state, but dint is 
smoother than Jiow, and this explains why the two cases pro¬ 
duce the same IR-optical flux. Indeed, dint soon overcomes 
Jiow, in regions where IR and optical frequencies are still 
efficiently produced. 

In conclusion, the good fits to the source SEDs that the 
inhomogeneous helical jet model can produce suggest that 
instabilities developing inside the jet and changing the jet 
geometry may explain the variations of the SED shape even 
when the intrinsic jet emission does not change. We notice 
that three-dimensional numerical simulations of relativistic 
magnetized jets show that kink instabilities are produced, 
leading to jet wi ggling and beam defl ection off the main 
longitudi nal axis jMignone et al 201Clli or he lical jet defor¬ 
mations (IMoII. Spruit fc Obergaulingerjl200^ . These results 
indicate that variations of the viewing angle of the jet emit¬ 
ting regions are likely to occur. 


6.2 Other possible explanations 

Another possible explanation for the optical-X-ray spectral 
shape and its variability may be the presence of an extra 
contribution to the absorption, in addition to the Galactic 


The laws according to which and d ecrease along the 
jet (see Eqs. 4 and 5 in IVillata &: Raiterilll999h are defined by: 
logZi = —2.3, logZ 2 = —1.8, Cl = C 2 = 3.3. The trend of emissiv¬ 
ity has been slightly modified with respect to Eq. 9 in Villata & 
Raiteri through the introduction of a further multiplicative factor 
(I/Iq)Oo/<^ 3 ) ^ where loglo = —1, cq = 3, C 3 = 2.3. 
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Figure 11. Spectral energy distributions of PG 1553+113 pre- 
sented in Figs. and 1101 fitted with the helical jet model by 
IVillata fc RaiteriTlIlQQd l. Three brightness states are presented: 
the high state of 2012 April 29 (dashed line), the intermediate 
state of 2013 May/July (dotted line) and the low state of 2009 
September/October (dot-dashed line). The fits were obtained by 
varying only the orientation of the jet, leaving the intrinsic flux 
unchanged. 
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Figure 12. The helical jet model used to fit the SEDs of Fig. 
111! Top: the trend of the maximum and minimum emitted fre¬ 
quencies as a function of the distance z along the helical jet axis. 
The behaviour of the emissivity is also shown (arbitrary units). 
Bottom: the viewing angle 6 and Doppler factor <5 along the jet 
axis for the three fits to the source SEDs presented in Fig. 1111 


one, which would imply less steep UV spectra and softer 
X-ray spectra, eliminating the need for an X-ray spectral 
curvature. This is questioned by the XMM-Newton data, 
which already present too high UV fluxes and for which the 
fit of the EPIC spectrum with a power-law model with Nn 
free is statistically worse than the fits with curved models. 

An alternative picture is that there are multiple (at least 
two) emission components, so that the optical-UV and X- 
ray emissions have a different origin. This would also ex¬ 
plain in particular why the PG 1553+113 X-ray flux varies 
by about a factor of 4 while the optical-UV level does not 
change, as can be derived when comparing the state of the 
source in 2009 September-October with that during the sec¬ 
ond MAGIC pointing (see Fig. [Toll . 

Another possibility is that the energy distribution of 
the relativistic electrons in the source, N{'y), where 7 is the 
Lorentz factor of the electrons, has a more complex shape 
than the usually assumed power-law or log-parabolic one. 

To distinguish among the various interpretations, we 
need more high-quality simultaneous observations in the 
optical-UV and X-rays, covering different brightness states 
of the source. This would allow us to better understand the 
correlation between the emission in the two bands. 


7 SUMMARY AND CONCLUSIONS 

We have presented multiwavelength data acquired during 
the WEST campaign organized in 2013 to study the blazar 
PG 1553+113 during MAGIC very high energy observations. 
We have complemented the radio-to-optical data from the 
WEBT with UV and X-ray observations by Swift and XMM- 
Newton to analyse the synchrotron emission in detail. A 
forthcoming paper will address the higher-energy emission 
and in particular the results obtained by MAGIC and Fermi 
(Ahnen et al. 2015, in preparation). 

The source revealed an enigmatic synchrotron be¬ 
haviour, the connection between the UV and X-ray spectra 
in the SED appearing more complex than usually assumed. 
While Swift-XK! observations have too low statistic to de¬ 
rive detailed information on the X-ray spectral shape, the 
long exposure by XMM-Newton in July 2013 confirms that 
the X-ray spectrum is curved, but probably not as curved 
as a log-parabolic fit would imply. On the other side, Swift- 
UVOT as well as HST-COS and lUE data show a steep UV 
spectrum that hardly connects to the X-ray one. 

Among the possible interpretations, we have investi¬ 
gated a scenario where the observed flux variations are due 
to changes of the viewing angle of the emitting regions in 
a helical jet, while the intrinsic flux remains constant. We 
have obtained good fits to the SEDs built with contempo¬ 
raneous data in different brightness states by varying only 
three parameters defining the jet orientation. We are not 
claiming that all the source variability is due to geometrical 
effects, but our analysis reveals how important orientation 
effects are in these beamed sources. 

We finally note that a periodicity of 720 ± 60 d has 
recently be en found when ana lysing the 7-ray data of PG 
1553+113 dGiprini et al.ll^lJh . The helical jet model can 
in principle explain periodicities by assuming that the he¬ 
lix rotates b ecause of the orbital motion in a binary black 
hole system dVillata et al.l[l9^ : lOstorero, Villata fc Raiteril 
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l2004h . However, a reproduction of both the SED and multi¬ 
wavelength flux variations requires good data sampling and 
a fine-tuning of the model parameters. 

A further observing effort needs to be spent to fully un¬ 
derstand the emission variability mechanisms of this object. 
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